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Abstract.  In  this  paper  we  consider  the  ignition  and  structure  of  a  reacting  compressible  mixing  layer 
using  finite  rate  chemistry  lying  between  two  streams  of  reactants  with  different  freestream  speeds  and  tem¬ 
peratures.  Numerical  integration  of  the  governing  equations  show  that  the  structure  of  the  reacting  How  can 
be  quite  complicated  depending  on  the  magnitude  of  the  Zcldovich  number.  In  particular,  for  sufficiently 
large  Zcldovich  number,  the  three  regimes  first  described  by  Linan  and  Crespo  (1976);  ie.,  ignition, 
deflagration,  and  diffusion  flame,  recur  in  supersonic  as  well  as  in  subsonic  flows.  An  analysis  of  both  the 
ignition  and  diffusion  flame  regimes  is  presented  using  a  combination  of  large  Zcldovich  number  asymptot¬ 
ics  and  numerics.  This  allows  us  to  analyze  the  behavior  of  these  regimes  as  a  function  of  the  parameters  of 
the  problem.  For  the  ignition  regime,  a  well  defined  ignition  point  will  always  exist  provided  the  adiabatic 
flame  temperature  is  greater  than  either  freestream  temperature.  One  important  result  is  that  at  supersonic 
speeds  ignition  occurs  far  downstream  from  the  plate  and,  as  the  flow  is  accelerated  to  hypersonic  speeds, 
ignition  is  exponentially  delayed.  For  the  diffusion  flame  regime,  the  location  of  the  flame  changes 
significantly  with  changes  in  the  equivalence  ratio  and  the  Schmidt  numbers. 


This  work  was  supported,  in  part,  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract  No.  NAS1-I8605  while 
the  authors  were  in  residence  at  the  Institute  for  Computer  Applications  in  Science  and  hnginccring,  NASA  I^anglcy  Research  Center, 
Hampton,  VA  23665. 
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1.  Introduction.  Quite  recently  it  has  been  realized  that  an  understanding  of  the  structure  and  the  sta¬ 
bility  characteristics  of  compressible,  particularly  supersonic,  reacting  mixing  layers  is  extremely  important 
in  the  context  of  the  scramjet  engine.  We  have  begun  a  systematic  study  (Jackson  and  Grosch,  1990)  of  die 
inviscid  spatial  stability  problem  of  these  Hows  with  the  combustion  zone  modeled  by  the  llame  sheet,  as 
f  described  by  Jackson  and  Hussaini  (1988).  We  found  that  the  addition  of  combustion  in  the  form  of  a  llame 

I  sheet  had  important,  and  complex,  effects  on  the  How  stability.  In  addition,  we  also  showed  that  with 

'  sufficient  heat  in  the  llame  sheet  the  flow  can  undergo  a  transition  from  convective  to  absolute  instability, 

even  without  the  presence  of  a  region  of  reversed  llow  in  the  layer.  There  arc  of  course  situations  in  which 
the  llame  sheet  model  may  not  be  a  sufficiently  accurate  representation  of  the  combustion  process.  This 
could  be  important  in  the  calculation  of  the  stability  characteristics  of  the  reacting  How.  Therefore  in  diis 
paper  we  return  to  the  problem  of  determining  the  structure  of  Hie  laminar  diffusion  flame  in  a  compressi¬ 
ble,  particularly  supersonic,  mixing  layer  in  using  finite  rate  chemistry. 

In  a  seminal  paper,  Linan  and  Crespo  (1976)  examined  the  structure  of  a  diffusion  llame  in  die 
i 

i  unsteady  mixing  of  two  half  spaces  of  fuel  and  oxidizer.  By  using  a  combination  of  large  activation  energy 

asymptotics  and  numerics  they  showed  that  three  laminar  regimes  exist  within  the  How,  the  ignition, 
deflagration,  and  diffusion  regimes.  These  same  regimes  were  also  found  by  Jackson  and  Hussaini  (1988)  in 
their  study  of  combustion  occuring  between  two  parallel  supersonic  streams  with  nearly  the  same  speed  and 
temperature.  These  three  regimes  arc  shown  in  Figure  1.  The  ignition  regime  is  a  region  where  die  com¬ 
bustible  gases  mix  until,  at  some  tinite  distance  downstream  of  the  plate,  a  thermal  explosion  occurs  and  die 
gas  is  ignited.  The  second  regime  is  the  deflagration  region.  After  ignition,  a  pair  of  well -defined 
deflagration  waves  (or  "premixed  tlameleis")  emerge  according  to  classical  thermal  explosion  theory  .  These 
waves  penetrate  Hie  mixing  layer  until  all  of  the  deficient  reactant  is  consumed.  Just  downstream  of  the 
deflagration  wave,  a  diffusion  flame  regime  exists  where  the  mixing  process  is  governed  by  diffusion  in  die 
direction  normal  to  the  flame.  This  regime  covers  the  bulk  of  the  overall  flow  field.  Collectively,  the  igni¬ 
tion  and  dellagrauon  regimes  can  occupy  a  very  small  portion  of  the  flow  field,  especially  at  supersonic 
speeds.  These  regimes  are  highly  nonlinear  and  hence  are  not  easily  amenable  to  analysis.  If  the  incoming 
flow  were  turbulent  raflier  than  laminar,  we  surmise  there  would  be  a  collection  of  initially  discrete 
premixed  llamcleLs.  Full  numerical  solutions,  based  on  the  compressible  Navier  Stokes  equations,  appear  to 
miss  Hie  ignition  and  the  deflagration  regimes,  and  hence.  Hie  premixed  Hamelet  is  not  captured.  Perhaps 
careful  fine-tuning  of  the  numerics  will  uncover  fliese  regimes. 

The  spatial  development  of  the  reacting  compressible  mixing  layer  is  considered  here  without  making 
the  assumption  of  infinite  Damkohler  number.  If  the  Danikohler  number  were  to  be  taken  to  be  infinite,  die 
flow  and  combustion  fields  could  be  described  by  similarity  solutions.  With  a  finite  Damkohler  number,  the 
flow  and  combustion  fields  evolve  in  Hie  downstream  direction  and  only  approach  the  similarity  solutions  of 
the  flame  sheet  at  infinite  distance  from  the  plate.  The  mean  fields  with  finite  Damkohler  number  will,  in  a 
future  paper,  be  used  in  stability  calculations. 

In  Hie  next  section  we  formulate  the  problem.  Numerical  solutions  to  the  full  equations  are  then 
given  in  secdon  3.  Sections  4  and  5  present  an  analysis  of  the  ignition  and  diffusion  flame  regimes  using  a 
combination  of  large  activation  energy  (and  hence  Zeldovich  number)  asymptotics  and  numerics.  Conclu¬ 
sions  are  then  given  in  the  final  section. 
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2.  Formulation  of  the  Mean  Flow  Equations.  The  nondimcnsional  equations  governing  the  steady 
two  dimensional  (low  of  a  compressible,  reacting  mixing  layer  with  zero  pressure  gradient  lying  between 
streams  of  reactants  with  different  speeds  and  temperatures  (Figure  1)  arc  given  by  (Williams,  1985) 


(pLO.  +(pF),  =0,  l=p7. 

(la« 

P  (UU,  +VUy)  =  {nUy)„ 

(lc) 

p(UTt 

+  I7y)  =  Pr-'diT,),  +  (y-  l),W2p^2+  pfh 

(Id) 

pm,.x 

+  V'Fy-y)  =  Sc,~l  (nFjwy)y  -  p,  Q.  j  =  1,2, 

(le) 

£2  =  D  pF,  F2e~y‘  r  . 

00 

In  these  equations,  the  x  axis  is  along  the  direction  of  (low;  the  y  axis  is  normal  to  the  flow;  U  and  V  are 

the  velocity  components  in  the  x  and  y  directions,  respectively;  p  is  the  density;  T  is  the  temperature;  and 

F  t  and  F:  are  the  mass  fractions,  with  the  reaction  assumed  to  be  irreversible  and  of  the  Arrhenius  type. 
The  viscosity  p  is  assumed  to  be  a  function  of  temperature.  The  nondimcnsional  parameters  appearing 
above  are  the  Prandtl  number  Pr ,  the  Schmidt  number  Sc,  =  Pr  Lc ,  for  species  j  where  Le,  is  the  Lewis 
number  for  species  j,  the  Mach  number  M  =  Ujam  ,  the  Zeldovich  number  Z c  =  E/R  7'„  with  E  the 

dimensional  activation  energy  and  R  the  universal  gas  constant,  D  the  Damkohicr  number  defined  as  the 

ratio  ol  the  characteristic  diffusion  time  scale  to  the  characteristic  chemical  time  scale,  p  the  heat  release 
per  unit  mass  fraction  of  the  reactant,  (3,  a  parameter  involving  stoichiometry,  and  finally  y  the  specilic- 
heats  ratio.  The  equations  were  nondimensionalized  by  the  freestream  values  /  „,  p,,,  (/„,  Flao  for  the  tem¬ 
perature,  density,  velocities  and  mass  fractions,  respectively.  Lengths  are  referred  to  some  characteristic 
length  scale  of  the  llow. 

Tiic  boundary  conditions  consistent  with  (1)  arc 

T  =  V  -  /'  i  =  1,  F 2  =  0  ai  x  =  0,  y  >  0,  and  x  >  0,  y  — >  (2a) 

f7 

I  -  P/ ,  U  =  p(,  <  I,  /•' ,  =  0,  F2  =  —  s  (p2 /  3i) 0-1  al  x  -  0,  y  <  0,  and  x  >  0,  y  — >  -«>,  (2b) 

'  i  ~ 

where  0  is  the  equivalence  ratio  defined  as  the  ratio  of  the  mass  fraction  F,  in  the  fast  stream  to  the  mass 
fraction  F2  in  the  slow  stream  divided  bv  the  ratio  of  their  molecular  weights  times  their  stoichiometric 
coefficients.  If  0  =  1,  the  mixture  is  said  to  be  stoichiometric,  if  0  >  1  it  is  Fj  rich,  and  if  0  <  1,  it  is  F, 
lean.  If  pr  is  less  than  one,  the  slow  gas  is  relatively  cold  compared  to  the  fast  stream,  and  if  pr  is  greater 
than  one  it  is  relatively  hot.  For  what  follows  we  take  the  stoichiometric  parameter  to  be  unity. 

The  mean  llow  equations  (1)  arc  first  transformed  into  the  incompressible  form  by  means  of  the 
Howarth-Dorodnitzyn  transformation 
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Y=jpdy,  V  =  pV  +  U\  p,  dy. 
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Next  we  transform  to 


T|  -  rj0 


Y 

2  \'I 


(4) 


which  is  the  similarity  variable  for  the  chemically  frozen  heat  conduction  problem,  and  r|„  correspond  to  a 
shift  in  the  origin.  For  the  case  of  both  streams  being  supersonic,  T)„  is  determined  uniquely  by  a  compati¬ 
bility  condition  found  by  matching  the  pressure  across  the  mixing  layer,  resulting  in 


ffp- 


V'(r,°°)  = 


-1 _ Pr_ 

2  P(2>  Pf1  -  1  Pr 


V  '(r  ,-o°). 


which  reduces  to  that  given  by  Ting  (1959)  for  =  1.  We  note  here  that  if  the  freestream  velocities  were 
subsonic,  then  the  compatibility  condition  would  be  trivially  satisfied,  and  thus  T|„  would  remain  indeter¬ 
minate  (Klemp  and  Acrivos,  1972).  Under  these  transformations,  with  U  -  f'W  and  U  -  (t|  f  f)  -x  , 
equation  (1)  become 


/"  + 2 //"  =  (),  (5a) 

4  JC  /  /  x  -  Pr~l  T"  -  Iff  -  (Y  -  1)  M\f")2  =  Ax  p£>  F^F2eVt  7  ,  (5b) 

4  xf’F;,x  -  Si.f'Fj  -  2  fFj  —  -Ax  D  F)F2c  ^  1  .  (5c) 

where  the  primes  indicates  partial  differentiation  with  respect  to  die  similarity  variable  T),  and  we  have 
assumed  the  linear  viscosity  law  p  =  T  (for  Chapman’s  linear  law'  p  =  C  7  ,  the  constant  C  can  be  scaled 
out  by  rescaling  T|  and  /  appropriately,  but  must  be  borne  in  mind  when  transforming  the  variables  back  to 
their  dimensional  forms).  In  terms  of  the  transformed  variables,  the  boundary  conditions  are 

l  =  f  -  F)  =  1,  F  i  =  0  at  x  =  0,  T|  >  0,  and  x  >  0,  r|  — >  «>,  (6a) 

F  =  P r<f  =  P(/ ,  F  |  =  0,  F 2  =  0_l  at  x  =  0,  T)  <  0,  and  x  >  0,  t|  — »  -°o.  (6b) 

Finally,  we  note  that  the  Damkohlcr  number  D  can  be  scaled  out  of  the  equations  by  a  rescaling  of  the  x 
coordinate.  Thus,  the  Damkohlcr  number  is  not  a  true  parameter  of  the  above  system.  However,  rescaling  x 
by  the  Damkohlcr  number  is  not  particularly  useful.  Since  D  is  typically  exponentially  large  the  rescaled 
coordinate  will  be  exponentially  stretched,  and  this  is  not  desirable  from  a  numerical  viewpoint. 

3.  Numerical  Solutions.  The  system  (5)  with  boundary  conditions  (6)  must  be  solved  numerically. 
Noting  that  /  is  only  a  function  of  T|,  wc  sec  that  (5b,c)  arc  parabolic  in  x  and  thus  can  be  solved  by  a 
marching  procedure  in  x.  Wc  do  this  by  a  Crank-Nicolson  procedure.  This  leads  to  an  implicit  nonlinear 
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system  of  equations  in  q  which  must  be  solved  at  each  x  step.  These  arc  solved  by  an  iterative  scheme 
using  the  boundary  conditions  in  q  as  given  in  (6). 

In  Figure  2a  we  plot  the  maximum  temperature  in  q  versus  x  for  Ze  -  10,  M  =  3t'  =  0, 

pr  =  p  =  q>  =  l.  The  adiabatic  flame  temperature  is  1.5  for  these  conditions  (sec  equation  (50)),  and  is  larger 

than  either  frccstrcam  temperature.  There  is  a  smooth  and  gradual  transition  from  the  inert  solution  at 
x  -  0  to  a  diffusion  flame.  In  Figure  2b  we  show  similar  results  but  with  Ze  =  30  and  all  other  parameters 
the  same.  Note  the  sudden  transition  between  the  inert  solution  and  the  diffusion  llame  solution  occurring 
in  a  narrow  range  in  x  around  x  =  3.  This  initial  region  defines  the  ignition  regime,  and  hence  the  liame 
location.  Similar  results  arc  shown  in  Figure  3  but  for  M  =2  with  an  adiabatic  flame  temperature  of  1.7. 
From  this  figure  it  is  apparent  that  the  ignition  regime  also  exists  in  supers onic  Hows.  As  the  Mach  number 
increases  there  is  a  corresponding  increase  in  the  inert  temperature  at  x  =  0  due  to  viscous  heating.  Ulus 
lessening  the  relative  effect  of  combustion  on  the  overall  temperature  field. 

In  Figure  4  we  plot  the  temperature  and  mass  fraction  profiles  as  a  function  of  q  at  various  x  loca¬ 
tions  for  the  same  conditions  as  in  Figure  2b.  In  Figure  4a  the  rapid  rise  in  the  temperature  over  a  narrow- 

range  of  x  should  be  noted;  this  indicates  ignition.  Note  also  that  the  temperature  peak  shifts  and  that  the 
profile  is  asymmetric  due  to  the  asymmetry  in  the  Lock  velocity  profile.  Figure  4b  shows  die  profile  ot  the 
F  ■  mass  fraction.  At  the  smaller  values  of  x,  F,  decreases  gradually  from  1  at  q  =  +«  to  0  at  q  - 
with  F •  varying  most  rapidly  over  the  range  -3  <  q  <  I.  This  shows  that  there  is  a  diffusion  of  F]  from 
the  q  >  0  region  into  the  q  <  0  region  w  ith  F,  =  0  only  at  q  -  -  ».  As  x  is  increased  the  rale  of  change  of 
F,  with  q  over  this  /one  increases.  For  larger  x  one  can  see  diat  F,  drops  to  zero  for  q  just  below  zero 
showing  the  presence  of  a  diffusion  llame  in  which  all  of  F,  is  consumed.  But  note  dial  there  is  a  small 
secondary  maximum  in  the  F.  distribution  in  —3  <  q  <  - 1 .  This  shows  the  presence  of  a  "premixed 
llame  let"  in  this  premixed  region.  The  distribution  of  F2  in  q  at  various  x  is  shown  in  Figure  4c  and  is  the 
converse  of  that  of  F,.  The  only  exception  is  in  the  secondary  maxima,  showing  that  the  presence  of  die 
"premixed  llame  let”  is  much  more  pronounced.  Again  the  asymmetry  of  the  secondary  maxima  is  due  to  die 
asymmetry  of  the  velocity  profile. 

The  existence  of  the  "premixed  llamclets"  and  the  diffusion  llame  is  shown  quite  clearly  in  Figure  5a, 
where  we  plot  the  loci  in  the  (x,q)  plane  of  the  maxima  of  £2,  the  chemical  production  term,  also  for  die 
data  of  Figure  2b.  As  shown  in  this  figure,  the  position  of  the  maximum  decreases  from  about  q  =  -0.2  at 
x  =  0  to  nearly  q  =  -0.6  at  x  ~  2.9.  The  maximum  value  of  £2  increases  along  diis  curve.  At  x  ~  2.9 
ignition  occurs  and  two  maxima  appear  giving  rise  to  the  "premixed  fiamelels".  Beyond  the  ignition  point 
the  premixed  flamclcts  move  outwards  in  the  shear  layer  until  all  of  the  deficient  reactant  is  consumed.  The 
appearance  of  the  third  maxima  just  behind  the  ignition  point  marks  the  appearance  of  the  diffusion  llame. 
As  x  is  further  increased,  the  diffusion  (lame  becomes  dominant  and,  as  x  — »  °°,  the  diffusion  llame  thins 
and  approaches  a  flame  sheet  characterized  by  local  chemical  equilibrium  (Jackson  &  Grosch.  1990).  Note 
that  the  diffusion  llame  location  is  constant  in  q  as  x  increases. 

In  Figure  5b  we  show  similar  results  for  the  same  values  of  the  parameters  except  that  3/  =  0.5  and 
3  =  1-5.  This  case  corresponds  to  unequal  frccstrcam  temperatures.  The  adiabatic  flame  temperature  is  1.5 
(again  see  (50)),  and  is  larger  than  cither  frecslream  temperature.  Note  that  the  ignition  point  has  moved 
into  the  region  of  higher  frccstrcam  temperature,  and  the  location  of  the  diffusion  flame  is  unchanged  while 
that  of  the  premixed  (lames  has  changed.  In  Figure  5c  we  show  results  for  3r  =  0.5  and  3  =  0.4.  The 
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adiabatic  flame  temperature  in  this  case  is  0.95,  and  is  smaller  than  the  frccstrcam  temperature  at  +  °».  In 
contrast  to  the  two  previous  figures,  there  is  no  well  defined  ignition  point;  the  premixed  flame  merges 
smoothly  into  the  diffusion  flame  whose  location  is  unchanged.  In  addition,  the  premixed  flamclcts  arc 
absent.  Finally,  we  note  that  the  behaviour  described  in  Figure  5  also  occurs  for  Mach  numbers  greater  than 
zero. 


In  the  following  two  sections,  an  analysis  of  the  ignition  regime  and  a  brief  discussion  of  the  diffusion 
flame  regime  is  given  using  a  combination  of  large  activation  energy  (and  hence  Zeldovich  number)  asymp¬ 
totics  and  numerics. 

4.  Ignition  Regime.  To  investigate  the  above  behaviour  analytically  for  large  Zeldovich  numbers,  we 
take  Le,  =  1,  Pr  =  1,  resulting  in  Sc,  =  1.  This  allows  us  to  consider  linear  combinations  of  (5b)  and  (5c) 
which  eliminates  the  source  terms.  The  solutions  of  the  resulting  equations,  using  the  boundary  conditions 
(6),  are  given  by 

T  +  p Ft  =  pr  +  (1  -  Pr  +  P)^+  -  Pt/)2lF(l  -  V).  (7) 

T  +  pFj  =  pr  +  p4>-'  +(1  -  Pr  -p0-')'F+  1^1  M\\  -  p,  )2T(l  -  S'),  (X) 


where 


S'  = 


=  <! 


(f  -  pr )/(l  -  Pr- ) 


(1  +  erf  t\)  /  2 


0  <  B,-  <  I 

Pc  =  > 


(.9) 


Thus  F  |  and  F2  can  be  eliminated  from  the  energy  equation  (5b),  yielding 

4  xf'T,  -  7  -  2/  7  -  (y  -  1)M2(1  -  pr)2^')2 

=  I  Pr  +0  -  Pr  +  P)H'  +  1^-Af2(l  -prV'S'O  -  S')-  7] 

P  * 

x  [Pr  +  P <J>-‘  +  (1  -  pr  -  P4TI)S'  +  A/2(l  -  Pr  )2 S' ( 1  -  S')  -  T )e~y*  '  (10) 

for  the  temperature,  subject  to  the  boundary  conditions 

7=1  at  x  =  0,  T|  >  0,  and  x  >  0,  T)  — >  »,  (1  la) 

T  ~  Pr  at  x  =  0,  q  <  (),  and  x  >  0,  r|  ->  -°°.  (1  lb) 

In  the  absence  of  chemistry,  equation  (10)  reduces  to 


4  xf'T,  -  /"  -  If  T  -  (y  -  1)W2(1  -  Pr  )2  (S'  )2  =  0, 


02) 
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which  possesses  the  inert  similarity  solution 

T,  ^  T,(r\)  =  1  -  (1  -  MU  -'»')  +  1y~MH  1  -  M2^0  -  'VI  (13) 

which  is  also  a  solution  of  (10)  for  small  x.  As  x  increases,  more  of  the  combustible  mixes  until,  at  some 
finite  distance  downstream  of  the  plate,  a  thermal  explosion  occurs  characterized  by  significant  departure 
from  the  inert. 

To  analyze  the  ignition  process,  we  expand  about  the  inert  by  defining 


T,  =  Ze  (T  -  T,).  (14) 

and  take  the  asymptotic  limit  Ze  — »  °°.  Substituting  into  (10),  we  see  that  the  left  hand  side  becomes 
()(Ze~])  while  the  right  hand  side  is  0(D  e  7"  T’),  where  7’.  =  max(Ti).  In  order  to  balance  these  terms, 
D  must  therefore  have  the  form 


D  = 


0  Z»  T. 

pSZc  e 


(15) 


where  5  will  be  chosen  in  the  course  of  the  analysis.  Substituting  (14)  and  (15)  into  (10)  yields  the  follow¬ 
ing  equation  for  /', 


4x  7  ,.,  -  /■:  -2/  7  ;  =4*-^ipV 

riP' 


Ze'1/, 


[P0-‘(1  -  'V  )  -  Ze-'Ti  I  exp 


[  7  ,  +  Ze  (7,  -  7  . ) 

[  T.  (7/  +  Ze'-1  /  ,)  j 


(16) 


which  must  be  solved  subject  to  the  boundary  conditions 

7  ,  =  0  ut  x  =  0  and  at  t)  ->  ±°°.  (17) 

Note  that  the  right  hand  side  of  (16)  is  exponentially  small  except  where  |7'/  -  T.  \  =  0( Ze"1).  There  arc 
now  several  cases  to  consider,  depending  upon  the  magnitude  of  the  parameters  Pj-,  P(/,  and  M .  Each  case 
is  described  below. 

4.1.  |1-P,|<<1,  |l-p„|<<l,  and  M  =  0(1). 

This  case  corresponds  to  the  simple  but  important  case  of  ignition  for  nearly  equal  frccstrcam  temperatures 
and  velocities,  and  the  analysis  is  presented  by  Jackson  and  Hussaini  (1988).  For  this  reason,  we  only  briefly 
review  the  results  here.  We  begin  by  setting 


Pr  -  Ze  (I  -  Pr),  P;/  -  VZe  (1  -  p(/). 


(18) 
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with  j3r,  and  M  fixed  and  0{  1).  In  particular,  we  see  that 


/  =  ti  +  0(Ze'V2),  V  =  V0  +  0(Ze-''2),  ¥0  =  -  (1  +  erfr\). 


(19) 


The  inert  solution  (13)  is  expanded  as 


T,  =  1  +Ze-‘[-pT(l  -1*.)+  4».(1  -  4 'a)l  +  0(Ze~2). 


(20) 


which  yields  7”.  =  1.  Substituting  into  (16),  together  with  the  choice  8=1,  yields  the  following  equation 
for  the  disturbance  temperature 


4 xT,  -  2  T)  7',  -  T,  =  4  jc  4>0  ( 1  -4%)  exp 


T1-It(1-'V0)+1y-M2U'VA1-'Vo)  ■  (21) 


subject  to  the  boundary  conditions  (17).  This  equation  was  solved  numerically  by  a  Crank-Nicolson  implicit 
scheme.  As  x  increases,  the  solution  for  T !  becomes  unbounded  at  some  finite  location  (x.  ,q. )  downstream 
of  the  plate,  indicating  thermal  runaway.  The  results  arc  presented  in  Figure  6,  where  we  plot  x.  and  q.  as 
a  function  of  M  for  various  values  of  $r.  Note  that  for  a  fixed  3r,  increasing  M  [5(I  causes  die  ignition 
point  to  move  closer  to  the  origin.  For  a  fixed  M  M  increasing  (or  decreasing)  j^r  from  zero  (i.c.,  cooling 
(or  heating)  the  gas  from  below)  causes  q.  to  move  into  the  region  of  higher  frccstrcam  temperature.  In 
addition,  x .  increases  as  the  gas  is  cooled  below  and  decreases  as  the  gas  is  heated  below. 

As  noted  above,  as  (5r  increases  or  decreases  from  zero,  the  ignition  point  moves  into  die  region  of 
higher  freestream  temperature.  This  case  corresponds  to  different  frccstrcam  temperatures,  and  can  be 
analyzed  in  the  limit  [M  — >  We  shall  assume  (37  <  1  (j>r  >  0)  and  hence  the  ignition  point  is  located  in 
the  region  q  >>  l.  The  analysis  for  the  case  >  1  ($T  <  0)  is  similar.  Following  the  analysis  of  Linan  and 
Crespo  (1976),  we  set 


Ml -¥,)  =  *,  0  <  z  <  °°,  (22) 

and  using  the  asymptotic  expansion  of  H^,  we  sec  that 

q2  =  In ($r  / 2  Vrt)  -  Inq  -  Inz.  (23) 

The  reaction  zone  is  located  where  z  =  0(1),  and  to  a  first  approximation  wc  get 

rj?  =  ln({5r/2V7t)  —  Inq.,  (24) 

where  q.  defines  the  ignition  location.  Since  the  chemistry  is  confined  to  a  thin  region  about  q. ,  we  shift  to 
the  reaction  zone  coordinate 


1 

2q. 


In  z. 


q  =  q.  - 


(25) 
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Substitution  into  (16),  together  with  =  1,  and  neglecting  terms  of  0(T).~2),  yields  a  canonical  equation  for 
/',.  Repealing  the  analysis  for  Pr  >  1  yields  a  canonical  equation  for  'l\  of  exactly  the  same  form,  the 
differences  being  in  the  definition  of  the  parameters,  the  flame  now  being  located  in  the  slow  stream  where 
the  gas  is  initially  hotter,  and  T.  -  pT.  This  canonical  equation,  originally  given  by  Linan  and  Crespo 
(1976)  in  a  different  context,  is 


-*i '  z2V«  =  U  -  Ay)ev  1 

where 

Pr  <  1 

Pr  >  1  - 

This  equation  must  be  solved  subject  to  the  boundary  conditions 


''(P/-JK) 


07  <  1 


r 


X  I  (-fir1)’)  Pr  >  F 


j  Od  -  Pr)/P 
(Pr  -  I  V  P 


(26) 


(27) 


y(0)  ---(),  y2  (»)  =  ().  (28) 

The  first  boundary  condition  is  a  restatement  of  (17),  while  the  second  boundary  condition  comes  from 
matching  with  the  chemically  lro/en  region  q  -  0(1).  Note  that  x  is  now  a  parameter  of  the  equation. 

Linan  and  Crespo  have  shown  dial  for  A  <  l,  y  is  doubled  valued  for  x,  <  xu,  single  valued  at 
r;  t;,.  and  does  not  exist  for  t;  >  X|..  Near  xtw  we  note  that  Oy/dx.  becomes  large  showing  that  die  x 
partui I  derivative  can  no  longer  be  neglected  in  (16).  Therefore  x1<  may  be  considered  as  the  ignition  point. 
For  A  >  1,  y  is  single  valued  for  all  x.  indicating  that  there  is  no  well  defined  ignition  point.  In  this  case 
the  premised  flame  merges  smoothly  with  the  diffusion  llamc  and  there  arc  no  premised  flamclets.  This 
behaviour  is  shown  numerically  in  Figure  5c.  In  Figure  7  we  plot  die  ignition  point  ,t1>  versus  A.  As  A 
approaches  one  from  below,  the  ignition  point  moves  off  to  infinity,  while  for  A  >  1  there  is  no  well  defined 
ignition  point. 

4.2.  !  1  -  P  /  i  <  <  1 .  Pi  <  1,  and  M  <  <  1. 

This  case  corresponds  to  nearly  equal  freestreani  temperatures  in  a  shear  flow  at  small  Mach  numbers,  and 
we  therefore  set 


py  =  Ze  (1  -  p; ),  M  =  vZe  M , 


(29) 


with  p,  and  M  fixed  and  f>(l).  The  inert  solution  (13)  is  now  expanded  as 

T,  =  1  +/c",[-0f(l  -M»)+  A/2(l  -  p{,)2T(l  -  y)l  +  C(Ze-2), 


(30) 


again  yielding  /'.  =  1.  Here,  'F  is  given  by  (9a).  Subsututing  into  (16),  again  with  the  choice  5=1,  yields 
the  equation  for  the  disturbance  temperature 
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4x  f  T\  -2/  '/',  -  7',  =  4jc  4^(1 -y)  exp 


V,  -  3,  (1  -  40  + 


—  Af2(l 


3ty  )2  41  (1  -  4*) 


,  0\) 


subject  to  the  boundary  conditions  (17).  This  equation  was  again  solved  by  a  Crank-Nicolson  scheme.  As  x 
increases,  the  solution  for  T,  becomes  unbounded  at  some  finite  location  (x. . rj. )  downstream  of  the  plate, 
indicating  thermal  runaway.  Some  results  arc  presented  in  Figure  8,  where  we  plot  x.  and  T|.  as  a  function 
of  3;/  for  various  values  of  M  and  3r  =  0.  We  first  note  that  for  p(/  =  1,  the  ignition  point  is  independent 
of  the  Mach  number.  For  a  fixed  M ,  decreasing  3u  from  one  causes  the  ignition  point  to  move  back 
toward  the  plate  and  into  the  slower  moving  gas.  At  fixed  3t/  <  I,  an  increase  in  M  causes  the  ignition 
point  to  shift  toward  the  origin. 

Increasing  or  decreasing  3/  has  the  same  effect  on  the  ignition  location  as  in  the  previous  case.  To 
see  this,  we  first  note  that  the  asymptotic  expansion  of  die  Lock  profile  (9a)  for  large  r|  is  given  by 


f'  =  \ 


1  -  A  (q  +  B  )  1  exp  (-  (q  +  B  )2) 

3,  +  A  (rj  B  P  exp  (  -  3,  (rj  +  B  )2) 


T)  -4  t  CO 

T|  — >  - 


(.12 


where  A  and  B  arc  known  constants.  Following  the  analysis  of  the  previous  section,  we  again  obtain  die 
canonical  equation  (26)  for  the  temperature  perturbation,  with  x,  now'  given  by 


x  ‘ (3;  (1*  f  B  y)  3/  <  1 


x  '(  -  3/  in-  +  b  v)  3/  >  i- 


(■VI) 


and  A  as  before. 


4.3.  3r>0,  0  <  3,  <  1 ,  and  M  >  0. 

For  Mach  numbers  greater  dian  zero,  3<  <  I  and  any  value  of  3/  .  the  inert  solution  (13)  has  a  maximum  at 
I )  =  0;  i.c.. 


4>'  1  -  3, 


3i.  ro  -  24>) 


=  o. 


(34) 


II  4>  =  0,  then  the  location  of  the  reaction  zone  will  occur  in  the  outer  region  |t|  |  >>  I  where  die  gas  is 
initially  hotter.  Since  the  analysis  for  this  case  will  be  similar  to  that  of  the  previous  two  eases,  we  w  ill  not 
present  it  here.  If  the  second  term  of  (34)  is  zero,  then  the  reaction  zone  will  occur  in  the  mixing  layer,  and 
is  given  by  the  implicit  relation 


4\ 


1  +  _ 1  ~  3r _ 

2  +  (y-  1)A/2(1  -  3,;)2' 


(35) 
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The  corresponding  maximum  value  of  7}  is  given  by 

1  —  Pt  y  —  1  7 

i=I - —  +  -I - -  M2 


T.  =  max  (Tj)  =  1 


2  8  M  (1  ^  +  2(y-  1)M2(1  -  p(f)2' 


(36) 


Substituting  (34-36)  into  (16)  shows  that  the  right  hand  side  of  (16)  is  exponentially  small  except  in  a  nar¬ 
row  region  where  |7)  -  7.  j  =  0(Ze~l).  Outside  this  region  the  equation  is  essentially  frozen,  and  is  given 
by 


4-*  7  i.i  =  7  i  +  2/  T i ,  (37) 

valid  on  either  side  of  the  reaction  zone.  To  study  the  structure  of  the  reaction  zone  we  define  the  new 
coordinate 


T|  =  T).  +  8z, 


(38) 


where  8  is  the  thickness  of  the  reaction  zone  and  is  not  yet  specified.  The  temperature  and  velocity  are  also 
expanded  as 


/■;  =  7  ,.<*)  +  8/  ,,(x,r|)  +  •  •  •  .  =  4».  +  8:  + 


(39) 


-Substitution  into  (16)  yields 


=  -Ax  (.1 


4*.  )exp 


-  PeW. 


.  ,  ,  ] 


T.2 


where  we  have  chosen  §  Zc  !  *.  Note  that  (40)  can  be  integrated  once  to  yield 


<7  7  n 


dz 


-d(x)  j 


-=r  expl  -  b2z2 1  dz  +ct(x). 


where 


d(x)  =  Ax  4N  (1  -  SMcxpl  l\J  /  .21 


\n 
b  ’ 


(40) 


(41) 


(42a) 


X_LM2(1  „p(,)2(vp;)27.-2 


(42b) 


Evaluating  (41)  at  z  =  ±»  yields  the  matching  conditions  for  (37).  Specifically,  we  now  solve  (37)  subject 
to  the  boundary  and  jump  conditions 


l\(x,°°)  =  0, 
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(43) 


I\{x  ,T)"+  )  =  T,.(x). 


dl , 

— -(x,r[’+)  =  c,(x)  -  d(x), 

dr\ 


I  i(x ,  -  °°)  -  0.  7* , (x . r| -  )  =  7\,(x), 


dl\ 

-7— (x.TT  ) 
d  r) 


Ci(x), 


(44) 


with  the  initial  conditions 


7-1(0.T|)  =  0.  p  *  p. .  (45) 

Equation  (37),  valid  on  either  side  of  the  reaction  /one,  can  be  solved  numerically  by  marching  forward  in 
x  The  unknown  functions  Tim(x  )  and  r'i(x)  are  chosen  at  each  x  location  so  as  to  satisfy  the  boundary  con- 
diuons  at  r|  =  ±°°.  Thermal  runaway  will  occur  at  (x.  ,p.),  thus  defining  the  llame  location.  In  Figure  9  we 
show  7  ;>  and  c  is  a  function  of  x  for  p7  =  0.5,  pr  =  0,  and  M  =  2.  Note  that  as  x  approaches  x. ,  1 
becomes  unbounded  indicating  thermal  runaway.  In  Figure  10  we  plot  the  ignition  location  x.  and  p.  as  a 
function  of  the  Mach  number  for  p(  ~  0  and  various  values  of  p r.  As  the  Mach  number  becomes  large  x. 
moves  off  to  tx  w  hile  p.  tends  to  zero  for  any  p; .  This  shows  that  in  the  supersonic  region  ignition 
occurs  far  downstream  of  the  plate  and,  as  the  How  is  accelerated  to  hypersonic  speeds,  ignition  is  exponen¬ 
tials  delayed.  For  unequal  freestream  temperatures  (Py  z  1),  as  the  Mach  number  goes  to  zero  then  p. 
moves  out  ol  the  mixing  layer  into  the  region  of  higher  freestream  temperature.  When  this  occurs  this 
analysis  is  no  longer  appropriate,  and  one  must  consider  the  case  of  4*  =  0  in  (.74). 

5.  Diffusion  Flame  Rev  •  \.  Typically,  a  diffusion  flame  is  characterized  by  a  chemical  reaction  time 
that  is  much  smaller  than  a  characteristic  diffusion  time.  Chemical  reactions  then  occur  in  a  narrow  zone 
between  the  fuel  and  die  oxidizer,  where  the  concentrations  of  both  reactants  are  very  small.  Mathemati- 
cally.  the  assumption  of  very  last  chemical  reaction  rates  leads  to  the  limit  of  infinite  Damkohler  number 
which  reduces  the  diffusion  llame  to  a  llame  sheet  (i.e.  local  chemical  equilibrium).  This  assumption 
significantly  reduces  the  complexity  of  the  problem  since  it  eliminates  the  analysis  associated  with  die 
chemical  kinetics.  For  many  flows,  the  assumption  of  local  chemical  equilibrium  adequately  predicts  die 
location  and  the  shape  of  the  diffusion  flame  (Buckmastcr  and  Ludford,  1982;  Williams,  1985).  Since  die 
llame  sheet  model  has  a  simple  similarity  solution  and  hence  is  amenable  to  analysis,  we  give  its  structure 
below.  For  finite  values  of  the  Damkohler  number,  equations  (5b,c)  must  be  solved  directly  using  die 
numerical  scheme  mentioned  above. 

Assuming  unit  Lewis  and  Prandll  numbers  and  a  linear  viscosity  law,  the  profile  associated  with  the 
flame  sheet  mode!  is  given  by  (Jackson  and  Hussaim,  1988,  Jackson  and  Grosch,  1990) 

7,  =  I  -  (1  +  <tr')(l  -  40,  F 2  =  0,  (46) 

7  =  P r  +  pr1  +  (I  -  Pr  -  P4L')4^  +  M\\  -  P,;)2lF(l  -  40,  (47) 

valid  for  p  >  p^ ,  and 
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F,=  0,  F2  =  <tf' -  (1  +<|r1)4'. 


(48) 


1  =  P7  +  0  -  Pr  +  P)H-  +  -L_  M2(l  -  P£/)2H<(1  -  4»), 


(49) 


valid  for  q  <  q^ .  Here,  q7  gives  the  location  of  the  (lame  sheet  where  both  the  reactants  vanish,  and  T 
takes  the  adiabatic  flame  value 


Tf  =  Pr  +  (1  -  Pr  +  P)^  +  -  Pr)2^  (1  -  Vf).  (50) 

The  flame  location  is  given  by  the  implicit  relation 

1 "  TT»  <5,) 

which  is  independent  of  p£  and  M .  This  is  only  valid  for  the  case  of  a  linear  viscosity  law  where  die 
momentum  equation  decouples  from  the  energy  equation.  For  a  more  general  viscosity  law,  the  (lame  loca¬ 
tion  will  depend  on  P7  and  M  as  well.  When  using  (.4)  to  transform  back  to  the  physical  coordinates,  the 
diffusion  flame  location  in  physical  space  will  depend  on  p(/,  0,  pT  and  M ,  independently  of  the  viscosity- 
temperature  relation.  The  internal  structure  of  the  diffusion  llame,  which  appears  as  a  discontinuity  on  the  q 
scale,  is  described  in  in  Jackson  and  Hussaini  (,1‘WX). 

The  above  assumption  of  unit  Lewis  and  Prandtl  numbers  led  to  a  closed-form  analytical  solution. 
Since  the  Prandtl  number  is  about  0.72  and  the  Lewis  numbers  are  not  one.  the  Schmidt  numbers  will  in 
general  not  be  one.  To  detenu inc  die  influence  of  die  Schmidt  numbers  on  the  location  of  the  diffusion 


llame,  we  begin  by  considering  (5b,c)  in  die  limit  of  infinite  Damkohler  number 

/V  1  7  "  +  2  fi  +  (y  -  1)  M2(f")2  =  0,  (52) 

Scf'  F’  +  2  fFj  =  0,  (53) 

valid  on  either  side  of  the  llame  sheet.  The  appropriate  boundary  and  jump  conditions  are  given  by 

/  ==  1,  Ft  -  I,  F;  =  ()  as  T|  -»  <*>.  (54a) 

I  Pr  f  f  p  Sc'1  F'  |  =  (),  |7  )  =  0,  F ,  =  F2  =  0,  at  q  =  q;  ,  (54b) 

/  =  pT  Ft  =0,  F2  =  0'1  as  q  — >  -<*>,  (54c) 


where  (  I  denotes  the  jump  (  )  -  (  )  at  q7 .  Integrating  the  equations  (53)  for  the  mass  fractions  yields 

F,  =  1  -ejfyf'dq. 
n 


F  2  =  0, 


(55) 
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valid  for  r|  >  t\^  ,  and 


f,=0.  F2  =  r'  -  c2  j  (Vf- ldr\. 


valid  for  q  <  T\f .  Using  the  boundary  condition  F,  =  F2  =  0  at  yields 


J  &)C'dT\ 

"f 


’V 


j  QV’)Scidn 


(56) 


(57) 


The  location  of  the  flame  sheet  x\/ ,  where  both  reactants  vanish  simultaneously,  is  given  by  the  implicit 
relation 


'V  .  - 

<J>Sc2  J  ('v  fC2dr]  =  Sc ,  ('Vf  )S‘: 2  Sr'  J  {H,)Sc'dr[.  <  5H> 

i/ 

obtained  from  the  jump  condition  (54b)  on  the  derivatives.  In  Figure  1 1  we  plot  the  llamc  location  as  a 
function  of  .Sc :  for  various  cases.  In  Figure  11a  we  show  the  variation  with  Sc2-  As  Sc:  increases  the  flame 
location  moves  into  the  slower  moving  stream.  In  Figure  lib  we  show  the  variation  with  the  equivalence 
ratio.  As  4>  increases  the  flame  location  also  moves  into  the  slower  moving  stream.  Finally,  in  Figure  1  lc  we 
show  the  variation  with  [lr.  As  jir  increases  the  flame  moves  into  the  faster  stream. 

6.  Conclusions.  In  this  paper  we  have  considered  the  ignition  and  structure  of  a  reacting  compressible 
mixing  layer  using  finite  rate  chemistry  lying  between  two  streams  of  reactants  with  different  Ireestream 
speeds  and  temperatures.  Numerical  integration  of  the  governing  equations  show  that  the  structure  of  die 
reacting  flow  can  be  quite  complicated  depending  on  the  magnitude  of  the  Zcldovich  number.  In  particular, 
for  sufficiently  large  Zcldovich  number,  the  three  regimes  first  described  by  Linan  and  Crespo  (1976);  ie., 
ignition,  deflagration,  and  diffusion  llamc,  occur  in  supersonic  as  well  as  in  subsonic  flows.  Our  numerics 
pick  up  the  premixed  flamclets  as  clearly  shown  by  the  distribution  of  the  production  term.  Again,  depend¬ 
ing  on  the  magnitude  of  the  Zcldovich  number,  there  can  be  a  gradual  or  sudden  transition  from  an  inert 
solution  to  a  diffusion  flame. 

An  analysis  of  both  the  ignition  and  diffusion  flame  regimes  was  presented  using  a  combination  of 
large  Zcldovich  number  asymptotics  and  numerics.  This  allowed  us  to  analyze  the  behavior  of  these  regimes 
as  a  function  of  the  parameters  of  the  problem.  For  the  ignition  regime,  a  well  defined  ignition  point  will 
always  exist  provided  the  adiabatic  flame  temperature  is  greater  than  cither  frccstrcam  temperature.  One 
important  result  is  that  at  supersonic  speeds  ignition  occurs  far  downstream  from  the  plate  and,  as  the  flow 
is  accelerated  to  hypersonic  speeds,  ignition  is  exponentially  delayed.  For  the  diffusion  flame  regime,  the 
location  of  the  flame  changes  significantly  with  changes  in  the  equivalence  ratio  and  the  Schmidt  numbers. 

There  arc  substantial  changes  in  the  temperature  and  mass  fraction  distribution  as  the  mixing  layer 
evolves  downstream.  These  arc  expected  to  have  a  considerable  effect  on  the  stability  characteristics  of  the 
flow.  The  changes,  and  particularly  the  rate  of  change,  in  the  mean  profiles  were  shown  to  be  very  sensitive 
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to  the  value  of  the  Zcldovich  number.  Therefore  we  expect  that  the  stability  characteristics  will  also  be 
quite  sensitive  to  the  value  of  the  Zeldovich  number.  The  stability  calculations  are  underway  and  will  be 
reported  at  a  later  date. 
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Figure  1.  Schematic  showing  the  reacting  mixing  layer,  where  the  adiabatic  flame  temperature  is 
greater  than  7',  and  /  2. 


